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Abstract 

Active galactic nuclei (AGN) are believed to be the source of ultra high energy 
cosmic rays (UHECRs, E > 10^* eV). Particles are assumed to be accelerated in 
the accretion disk and the plasma jets, produced due to conservation of angular 
momentum, to the highest energies, where they interact with each other and 
produce pions, which decay among others in neutrinos. 

For a known cosmic ray spectral behavior, the main parameters in the calcula- 
tion of the neutrino flux from proton-proton interactions are the target density 
riH and the ratio of electrons to protons fc- Using most recent neutrino flux 
limits from IceCube point source searches, we set limits on the target densities 
for 33 FR-I galaxies. The densities are shown to be smaller than 30 cm~^ to 
2-10'^ cm~^, depending on the source and when using a fixed electron to proton 
ratio of fe = 0.1. This implies that some cosmic ray acceleration sites, espe- 
cially those close to the core of the AGN, can already be excluded, or else that 
the ratio of electrons to protons deviates significantly from the commonly used 
value of 0.1. 

For Centaurus A (Gen A) and Messier 87 (M 87) we use Fermi observations to 
model the 7-flux, the neutrino flux and the resulting target density. The detec- 
tion of these neutrinos will help to find information about acceleration processes 
in the source. 

Keywords: Active galactic nuclei, FR-I galaxies. Inelastic proton-proton 
interaction. Target density, Gentaurus A, Messier 87 



1. Introduction 

Active Galactic Nuclei (AGN) are the most powerful permanent objects 
known. The observed luminosities range from 10^^ erg s^^, measured for nearby 
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galaxies, up to lO'*^ erg s~'^ for distant galaxies. The prevalent picture is a su- 
permassive black hole (SMBH), located in the center of the host galaxy, with 
gravitational energy as the source of the luminosity. 

The emission is spread widely across the electromagnetic spectrum, often peak- 
ing in the ultra-violet, but with significant luminosity in the X-ray and infrared 
bands. The emitted power varies on time scales of years, days or minutes. Due 
to angular momentum conservation, plasma is ejected and forms bipolar colli- 
mated jets, strong radio sources if the host galaxy is elliptical, or weak radio 
sources if the host galaxy is a gas rich spiral. AGN are already observed and 
identified as sources of high energy 7-rays and additional observations indicate 
that the arrival direction of the highest energetic CR might correlate with the 
position of Cen A and M 87 P, i]. 

High energy neutrinos in AGN are produced in coincidence with high energy 
7-rays when high energetic protons interact with a local target and produce 
secondaries, which decay among others into neutrinos. The exact mechanism of 
the energy and momentum transfer in AGN is still under debate. Matter which 
is attracted by the SMBH, the central engine of the AGN cannot directly fall 
into the black hole, since it possesses angular momentum. A disk of matter is 
formed within which magnetic viscosity transfers angular momentum outward 
and mass inward. In the core region protons can be accelerated to high energies 
via reconnection of the magnetic field to high energies Q . The highly energetic 
protons can interact with other protons and photons, producing secondary par- 
ticles as high energetic neutrinos, which can leave their point of origin almost 
unimpeded to be detected on Earth. A detailed approach to proton-photon and 
proton-proton interaction in the jet can be found in Q . 

Here we focus on proton-proton (p-p) interaction in the core region and concen- 
trate on the muon neutrino flux calculation for 33 FR-I galaxies. Assuming that 
the electrons lose all their energy to synchrotron radiation, the radio luminosity 
is equal to the electron luminosity 



Furthermore protons and electrons are expected to be accelerated at the same 
site, meaning that the proton luminosity can be determined by assuming a 
constant ratio /e = Lc/Lp between radio and proton luminosity [5|. Using this 
assumption, the proton luminosity can be estimated from radio observations of 
individual sources. A further astrophysical parameter is the target density nn, 
which determines the optical depth Tpp of the p-p interaction. It is given by 



where /q is the initial intensity and / the observed one. For the calculation the 
optical depth is given by 
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where nn is the target density, crinei(-Bp) the inelastic cross section for proton 
proton interaction and R the is the size of the interaction region. For the consid- 
ered parameter space (emission from the core region of the AGN, R < 100 kpc, 
with densities riu < 2 • 10'^ cm~^), Equation ^ can be approximated by a linear 
behavior of the optical depth. 

For Gen A and M 87 there are detailed 7-ray observations available, which will 
be used to normalize the 7-flux and to determine the neutrino flux and the re- 
sulting target densities nn at a given fraction /e @. 

This paper is constructed in the following manner. In section 2 the main calcu- 
lation will be introduced. In section[3J the results for the neutrino fluxes and the 
derived limits on the target density for 33 FR-I galaxies will be presented. The 
observations of 7-rays from Gen A and M 87 provides the possibility of deriving 
an exact value for nu, which is done in section |4l Gonclusions are presented in 
section [S] 



2. Modeling neutrino and 7-ray spectra from FR-I using the sources 
radio luminosity 

High energy neutrinos can be produced through inelastic proton-proton in- 
teractions, where high energy protons interact with ambient protons and pro- 
duce pions which decay into neutrinos. The muon neutrino flux at Earth is 
given by 0] 



V'"{Ei,) = Cose • C riH / CTincl 



diV, 



da; 



^{E^/x)F^{x,E^/x)^ with x = E^/E^ 
p 



(4) 

Here, Eosc considers the oscillation, nn is the target density in cm~^, (linoi is the 
inelastic proton-proton cross section in mb, given by 



CTincl = (34.3 + 81.88L + 0.25L^) 



Eth 



with L = In 



1 TeV J ' 
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with E'th the threshold energy for the 7r+ production. Neutrinos at the sources 
are created with the ratio 



((i^c + t^o) + + + {vr + Vr)) = (1:2:0) 



(6) 



Due to oscillation of neutrinos from the distant source to Earth the ratio to be 
expected here is 



((t/e + V^) + [v^ + V^) + {vr + Vr)) = (1:1:1), 



(7) 



meaning eosc = 1/3 for our calculation Q. 

The function F^{x, E^/x) gives the muon neutrino spectrum for a fixed proton 
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energy Ep and is divided into three summands, Fi, F2 and F3. The first one 
denotes the muon neutrino spectrum of neutrinos from direct pion decay, F2 
gives the spectrum of neutrinos produced by muon decay and ^3 considers the 
electron neutrinos produced 

Further, dNp/{dAdEp), the incident proton spectrum in units of cm~^ TeV~^, 
is given by 

^4|oxpf-^V (8) 



dAdEp EP \ Eq 

Here, Ap is the normahzation of the spectrum in units of cm^^TeV"^, Eq is 
the cut-off energy and p is the spectral index. In the following calculations, 
we use p = 2. While the cosmic ray spectrum might indeed deviate from an 
E^^ behavior (see e.g . M , IceCube point source limits are usually given for an 
spectrum only [lO|. In future work, it will be interesting to investigate 
the effect of different spectral indices as well, when IceCube limits are provided 
for other cases as well. For the two sources Cen A and M 87, where 7-ray 
measurements indicate a deviation from an behavior, we actually do use 

the observed values, as we do not rely on IceCube limits in that case. The 
cosmic ray normalization is connected to the total cosmic ray energy Wp via 

-Ep dEp - —j-j-r^ (9) 



dAdEp " " indLiz) 




^ (10) 
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47rdL(z) 
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The cosmic ray emission is isotropic, meaning that only the fraction (47r(iL(z)^)~^, 
with dL representing the redshift dependent luminosity distance, reaches the 
Earth. Due to the scaling law of the cross section with the nuclei number, the 
results apply for different cosmic ray compositions. We perform the calculations 
for protons only for simplicity, but expect the same results for a heavier com- 
position. The minimum energy of the cosmic rays is the threshold energy for 
the 7r+ production, E"™'" « 1.2 • 10~'^ TeV, and E'^^^ is the maximum cosmic 
ray energy. E'^'"^ w 10^ TeV. 

The total proton energy of a single source is given by 



Wp^ J Lpdt = LpJ dt (13) 
1 r dz' 
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In Equation (1131) we assume a constant proton luminosity Lp. The parameter 
is the Hubble time and E{z) is given by [ll| 



E{z) = y/nM{l + + W + + ^A- (15) 

In Equation (|14p we assume a ratio /c between radio and proton himinosity, 
given by 



-t^radio J dB,, 

"div; 



(16) 



Considering the mentioned assumptions the normahzation is 

^Earth ^ ^radio ^ 1 ^ (I7) 



3. Limits on the target densities for 33 FR-I galaxies 

Fanaroff and Riley divided radio galaxies according to the correlation be- 



tween morphology and luminosity into two groups, FR-I and FR-II [I4I. FR-I 
are brightest towards the center, while FR-II are brightest at outermost part of 
the jet. Later it was discovered that this behavior is correlated with a critical 
radio luminosity at 178 MHz, Lns = 2 • 10^^ W/Hz, dividing the radio sources 
in FR-I and FR-II, where FR-I have lower and FR-II higher luminosities [l3| . 
Detailed observation showed that the morphology is correlated to the energy 
transport in the sources. FR-I sources have bright knots along the jets, while 
FR-II have faint jets but bright hot spots at the end of the lobes, indicating that 
they appear to be able to transport energy efficiently to the ends of the lobes. 
FR-I sources on the other hand are inefficient in the way that a large amount 
of energy is radiated. Considering the AGN unification scheme, FR-I galaxies 
are assumed to be the misaligned counterparts of BL Lacs [l^, meaning that 
the non thermal beamed emission from the relativistic jets should be present in 
radio galaxies. Observation show that the ratio of nuclear luminosities of FR-I 
and BL-Lacs show a correlation with the orientation of FR-I galaxies, support- 
ing the assumption that they are correlated with each other. 
In this paper we focus on FR-I galaxies presented [l3|. The selection includes 
33 FR-I galaxies with a given radio luminosity Lradio at 178 MHz. The redshift 
ranges between z = 0.0037 — 0.29. Using the radio luminosity, the redshift z, 
considering the cosmology {Hq — 75 km s~^Mpc^^ flu — 0.27, VIa — 0.73) 
and fe the normalization Ap and the muon neutrino fiux are calculated. The 
latter can then be used to calculate the neutrino flux in dependence of the target 
density nn. 

3.1. Results 

Table [T] provides the upper limits for the target densities of the 33 FR-I 
galaxies presented in [IS^. Here, the radio luminosity of each single source, the 
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cosmology dependent parameters, luminosity distance dj,{z) the lookback time, 
and the constant ratio /e = 0.1, were used to compute the normalization Ap. 
The resulting muon neutrino fluxes are calculated and matched in normalization 
to the IceCube limit, resulting in an upper limit to the target density. 
Including the knowledge on neutrino flux limits provided by IceCube [loj allows 
thus to set an upper limit on the target density. As IceCube limits are only 
provided directly for a fixed number of sources, we use the declination dependent 
sensitivity from [lo| . Since we have detailed information about the high energy 
part of CRs, we can use observations of data to calculate the proton target 
density. The two assumptions we use is that a) UHECRs originate from AGN 
and b) protons are accelerated in the jet in the same way as electrons. Taking 
this assumption we calculate tt-h, which ranges between 20 cm^'^ for 3C 028 and 
1500 cm-3 for 3C 386 for /c = 0.1. 

Comparing our results with the model of Kazanaz and Elliason [l^, where a 
spherically symmetric accretion shock, accelerating a fraction of the inflowing 
plasma to the highest energies is considered, the proton density is given by 

nn{x) « 1.15 ■ 10^-^a;-3/2 cm'^. (18) 

Here, m M/{1 Moyr"!) is the accretion rate, Mg := M/(1O^M0) the black 
hole mass in units of 10^ solar masses. The parameter x := r/rs, with the 
Schwarzschild radius rs, gives the radial distance. Comparing Equation (|18ll 
with our results indicates that accelerated protons might originate from a max- 
imum orbit Xmax ~ 3 • 10"* for m = 0.1, from Xmax ~ 5 ■ 10** for m — 0.2 and 
from Xmax ~ 9 ■ lO'' for rh — 0.5 and from a minimum orbit Xmin ~ 2000 for 
TO = 0.1, from Xmin ~ 3000 for to = 0.2 and from Xmin ~ 6000 for m — 0.5. 



Considering the jet disk model [16|, [17[, the target density in the jet in the 
observers frame, is given by 

nn = 11 • ri46 Qj/iZp^ cm"^. (19) 

Here, L^q is the disk luminosity in units of 10"*^ erg/s, qj/i is the ratio between 
jet power and disk luminosity, Zpc is the distance from the origin in parsec 
and r is the boost factor of the plasma. The parameters for this model are 
r ss 10, L46 < lO^'^ and qj/i w 0.15 for Zpc < 1. The resulting densities of 
riH < 1700 (Zpc/0.1)~^ cm~^ are well consistent with the limits derived here 
and would only start to become inconsistent for Zpc ^ 0.1, a distance smaller 
than 10^ Tg. So far, it is not clear, how well the jet disk scenario works for FR-I 
galaxies, though, as their accretion disks are extremely faint and diSicult to 



observe [18 1 



4. Cen A and M 87 

For Cen A and M 87 we use the Fermi LAT observations to normalize the 
7-ray flux, which fixes the neutrino fiux directly. Thus IceCube limits are not 
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Scatter plot of the target density tin 



10,0001 





Figure 1: scatter plot of nH/(/G/0.1) as a function of sin(RA) (top) and the muon 
neutrino flux for 3C 028 (bottom). The red (dashed) lines in the scatter plot give the 
density range. The red (dotted) line gives the sensitivity of IC40. We show the muon 
neutrino flux for one source. The other sources have the same shape for the flux but 
a different value for the sensitivity. 
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Name 


sin(RA) 


sensitivity X 10 ''^ 
in TeVcm~^s~^ 


Normalization in 
TeV~^ cm~^ 


nH/(/o/0.1) 

— 3 

m cm 


3C 028 


0.24 


3.34 


287.16 


<23 


3C 029 


0.25 


3.34 


101.18 


<65 


3C 031 


0.29 


3.54 


48.93 


<149 


3C 066B 


0.58 


5.27 


88.80 


<140 


3C 075 


0.70 


6.42 


40.52 


<420 


3C 76.1 


0.72 


6.82 


63.21 


<420 


3C 078 


0.73 


7.09 


91.76 


<280 


3C 083.1 


0.76 


7.68 


109.26 


<280 


3C 084 


0.77 


7.53 


112.96 


<280 


3C 089 


0.80 


8.15 


296.93 


<93 


3C 264 


0.07 


3.28 


77.80 


<93 


3C 270 


-0.08 


7.15 


58.71 


<187 


3C 272.1 


-0.11 


9.46 


9.76 


<1027 


3C 274 


-0.13 


13.03 


457.65 


<47 


3C 277.3 


-0.23 


26.23 


98.58 


<280 


3C 288 


-0.42 


52.85 


348.95 


<187 


3C 293 


-0.47 


68.55 


84.06 


<700 


3C 296 


-0.56 


90.74 


45.78 


<1400 


3C 305 


-0.67 


117.76 


95.14 


<747 




-0.72 


i.'iO.i!) 


licS.69 


<.iS7 


3C 314.1 


-0.74 


132.83 


144.95 


<747 


3C 315 


-0.75 


143.88 


224.66 


<373 


3C 317 


-0.76 


141.05 


241.42 


<397 


3C 338 


-0.92 


206.47 


216.36 


<467 


3C 346 


-0.95 


223.69 


175.10 


<607 


3C 348 


-0.96 


215.00 


5529.19 


<23 


3C 386 


-0.99 


257.34 


68.99 


<1552 


3C 424 


-0.74 


143.88 


200.83 


<467 


3C 433 


-0.63 


100.33 


699.21 


<117 


3C 438 


-0.52 


72.81 


896.32 


<93 


3C 442 


-0.44 


59.59 


70.93 


<968 


3C 449 


-0.38 


43.26 


38.05 


<1073 


3C 465 


-0.09 


7.60 


185.54 


<58 



Table 1: The observed 33 FR-I galaxies. Column 1 gives the the names, column 2 the sine 
of right ascension, column 3 gives the sensitivity of IC40, column 4 the normalization and 
column 5 the target density for /e = 0.1 



8 



needed in this case. The Large Area telescope is a pair conversion 7-ray tele- 
scope, covering the energy range from 20 MeV to more than 300 GeV [6| . Due 
to its vicinity to Earth, Cen A has been well studied over the entire electromag- 
netic spectrum, from radio to 7-rays. Observations performed by experiments 
like the Auger observatory indicate that the origin of the highest energy CRs 
{E > lO^^eV) could correlate with the angular position of Cen A (l9| . 
M 87 one of the nearest {d = 16 Mpc) and best studied radio galaxies, just 
like Cen A is known for its bright arcsec-scale jet. It contains an SMBH with a 
mass of 21(3 - 6) • 10^ M©. At TeV energies M 87 is detected by H.E.S.S. 
MAGIC [m and VERITAS [22^. 



Due to the short observation period of MAGIC, VERITAS and H.E.S.S. in com- 
parison to Fermi LAT we focus on Fermi observations for M87. A further reason 
is the time variability of the 7- flux for E > 730 GeV (2^ , while Fermi gives due 
to the longer observation period time averaged fluxes. A satisfactory answer to 
the question of the source of high energy 7-rays is not yet found. 
For Cen A we also focus on Fermi observations, since Fermi LAT can resolve 



the inner region of Cen A [23|. A further reason is the 77 absorption, making 
it unlikely that the 7-ray flux measured by H.E.S.S. and Fermi, originates from 
the same region [ij. 

First we give an estimate on the target density for the two sources, by assuming 
that 7-rays are produced in the inner region, leading to the condition 

A = = Dsource- (20) 

Here, A is the mean free path and Dgourcc is the diameter of the source. Using 
this assumption excludes the other AGN from [6] , since the origin of the observed 
7-rays cannot be absolutely restored to the inner region of the sources. 
For the observed energy range, the cross section remains almost constant see 
(Equation ^) 

A=— (21) 

(22) 

30 mb nn 

— (^)^ ,3) 

Assuming that the cores are spherical, using the distance to the source, the 
diameter of the cores can be estimated, as 

D Cen A = 1.5-105 PC (24) 

D M 87 = 2.3 • lO'^ pc. (25) 
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For M 87 the 95% confidence error radius r95% — 0°.086 is used [2^. Considering 
our assumptions, the mean free path has to fulfill the following relation, 



A = D 

source 

(26) 

^ (/c/0.1)"^ = 70 cm-3 for Cen A (27) 
^ riH ifjO.iy^ = 430 cm-3 for M 87. (28) 

It has to be noticed that the assumptions used here are simplifications. For 
example, the geometry can be very complex due to physical processes, like 
clumping, so that the assumption of spherical geometry might not be sufficient 
to describe the real shape of the source. This has an influence on the target 
density. Equation (1261) . Furthermore we assume 7- production only in the inner 
region. Nevertheless the above calculation provides a first order approximation 
of the required target densities. 

We compare the results from the estimate with the results, using our model 
explained in section 2. Therefore, we use the radio luminosity of the inner region 
to calculate the normalization. The fraction 1/10 = icorc°/^tot'a°' observed for 
Cen A given in [2^ is assumed to be the same for M 87. Using a power law 
function, with F for the spectral index and Eq for the cut-ofF energy, a power-law 
spectrum of primary cosmic rays is used to model the 7— spectrum: 

^^^ApE-^exp{~E/Eo) (29) 

Here, we used Eq = lO^TeV as the maximum energy, which is needed to explain 
the observed flux of ultra-high energy cosmic rays. The Fermi data for Cen A 
and M 87 require spectral indices of FconA — 2.5 and Fm87 — 2.21. Considering 
the cosmology, the radio core luminosity and the distances, we calculate the 
normalization ^p, the 7-ray flux (Figure [2]) and the densities for /o=0.1: 

^P^^ ^ = 459 TeV^^ cm'^ => tih (fjO.iy^ = 10 cm^^ (30) 

Al^ 87 = 136 TeV"^ cm'^ ^ nn (/o/O.l)"^ = 35 cm'^ (31) 

As can be seen for Cen A the density estimated from the core size. Equation 
((27|) and the one calculated within the model presented. Equation ((30|) are in 
agreement. 

For M 87 the situation is different. Uncertainties like the unknown core size 
or the unknown core luminosity have a large influence on the estimate. To 
eliminate these uncertainties, the core size has to be resolved and the fraction 
of the core to the whole radio luminosity has to be known. For this purpose, 
improved measurements for the core size of M 87 will help. Additionally, the 
detailed observations of a larger sample of AGN in the future can contribute to 
receive a statistical sample of the fraction of core to total luminosity of FR-I 
sources. If the variance is of physical nature, on the other hand, the ratio /e 
might deviate significantly from /o — 0.1. 
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Figure 2: Gamma ray flux for M 87 (top) and for Cen A (bottom). The (red) points give the 
Fermi observation Q. The spectral index for the proton spectrum is F. 
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4-0.1. The density and muon neutrino spectrum for Cen A and M 87 

The muon neutrino flux can be derived directly from the fit to the 7- 
measurement and compared with the limits given in [lo|. As one can see, the 
flux is below the current sensitivity of the IC40 for both sources. It has to 
be mentioned that the uncertainties in the assumptions are significant, so the 
results have to be interpreted carefully. For M 87 the sensitivity is about one 
order of magnitude and for Cen A two orders of magnitude higher than the flux 
in the corresponding energy range. Long-term observations with IceCube and 
possible future extensions might be able to yield a first significant signal from M 
87. In the case of Cen A, the measured flux is simply too steep to be obesrvable 
by high-energy neutrino telescopes. Other emission regions, like outer parts of 
Cen A, might be more interesting and existing H.E.S.S. measurements reveal 
the possibility of such sources with a rather flat spectrum [2^ . Future and more 
sophisticated experiments like KM3NET will help to obtain better observations, 
since it is optimized for the southern hemisphere and therefor sensitiv to lower 
TeV energies [27 1. 

Since the neutrino flux is derived from the 7-measurements, uncertainties in 7- 
observations have an influence on the estimates of the neutrino flux and hence 
an influence on the difference between calculated neutrino flux and the sensitiv- 
ity of the detector. The more precise the 7-observations are, the more precise 
are the calculated neutrino fluxes, resulting in an improved comparison between 
sensitivity and the flux. For this reason experiments with improved equipment 
like more sophisticated satellites and the future ground based Cherenkov de- 
tector array (CTA) will help to obtain more detailed 7-ray observations. If the 
spectrum of Cen A turns out to be signiflcantly flatter than it is currently indi- 
cated by Fermi measurements, the central source might still be interesting for 
neutrino telescopes. 
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Muon Neutrino Flux Of M87 




Figure 3: Muon neutrino flux for M 87 and Cen A. The black (solid) curve is the spectrum 
and the red (dashed) the electron neutrino spectrum. The dark blue (dotted) curve is the 
muon neutrino spectrum, if oscillations are considered @| . The blue (dashed) horizontal line 
gives the IceCube sensitivity for M87 at the northern sky and for Cen A at the southern sky 
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5. Conclusions 

In this paper, we derive limits on the product ?t.h//c/(0.1) for 33 FR-I galax- 
ies, using the radio luminosity at 178 MHz and neutrino flux limits. The target 
density can be limited to be between 20 cm^'^ and 1500 cm^'^, for a fixed /e = 0.1 
for the sources. An explanation is the normalization of the sources. Two ob- 
jects at the same location at the sky but different normalization, differ in the 
target density, since it is A^, oc nn, resulting in A^^/Ai,2 — n-ffa/^-Hi fo'" the two 
objects. Two objects with a different location but approximately same normal- 
ization, the ratio is nHi/nH2 c< cki/a2, with a for the sensitivity of the detector 
for the object. Considering different models of the distribution of matter in 
AGN, we conclude that the innermost core models for proton acceleration can 
be excluded, as densities there are expected to be of the order of 10^ cm~'^. An 
alternative explanation would be that the ratio of protons to electrons must be 
significantly lower than expected from standard theory 0, [2^ . 
Gamma-ray observations exist for Cen A and M 87, the two closest and best 
studied AGN and they are used in this paper to estimate the average density of 
the emission region in case of hadronic interaction processes. For our calculation 
we used Fermi observations, since the observation period by Fermi reflects the 
average flux over a long-time measurement, rather than individual short-term 
measurements by Imaging Air Cherenkov Telescopes like MAGIC, VERITAS 
and H.E.S.S. Those provide interesting insights on the flaring behavior of the 
sources, but no reliable measurement of an average flux at this point. We find 
that the density required for gamma-ray emission in the core of Cen A and M 
87 has to be of the order of 10 — 100 cm~^. The exact value again depends on 
the ratio of protons to electrons accelerated at the source. 

In the first part we assumed that 7-rays are only produced in the inner 
region and that the inner region is approximated by a spherical symmetry. We 
estimated the target density and compared these results with the densities we 
obtained by using the model explained in section 2. For Cen A the differences 
between the estimate and the model are in agreement within the uncertainty 
of the measurements, while for M 87 the situation is different. The differences 
can explained by our assumptions. Due to physical process, like clumping, the 
geometry can be very complex. Additionally 7-ray production might not only 
takes place in the inner region of the objects. In the case of M 87 neither radio 
luminositiy observations of the core region nor observations of the core region 
exist. This leads to a bigger difference between estimate and model. Future 
measurements will give more detailed information about the inner region of the 
sources and thus help with the identification of the sources of ultra-high energy 
cosmic rays. 
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